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We discuss a new method to determine the low-energy couplings of the AS — 1 weak Hamiltonian 
in the £ -regime. It relies on a matching of the topological poles in 1 /m 2 of three-point functions 
of two pseudoscalar densities and a four-fermion operator computed in lattice QCD, to the same 
observables in the Chiral Effective Theory. We present the results of a NLO computation in 
chiral perturbation theory of these correlation functions together with some preliminary numerical 
results. 
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1. Introduction 

Non-leptonic kaon decays can be described in the context of a Chiral Effective Theory (ChPT), 
where the non-perturbative dynamics is encoded in the low-energy couplings associated with the 
AS = 1 Hamiltonian. These low-energy couplings can be determined from first principles in lattice 
QCD, by performing a matching of suitably chosen correlation functions computed in lattice sim- 
ulations and in the effective theory [jjj]. Such matching should be carried out as close as possible to 
the chiral limit. 

In [Q], a new strategy to reveal the role of the charm quark mass in the AI = 1/2 rule was 
proposed. The idea is to consider the SU(4)-flavour limit, that is a theory with four light quarks 
corresponding to the u,d,s and c, and compute the low-energy couplings of the corresponding 
AS = 1 Chiral Hamiltonian. In a second stage, the charm quark mass is increased towards its 
physical value, monitoring the change of the low-energy couplings. In the first determination 
of the leading-order low-energy couplings of the SU(4) AS = 1 Hamiltonian, g , was presented. 
For details on this computation and the precise definition of these couplings we refer the reader to 
references [|||3|]. 

In this work we present a new method to determine g 1 * 1 in the so-called £ -regime [|], H], from 
correlation functions involving pseudoscalar densities that contain topological poles in l/(mV) n , 
when evaluated in sectors of fixed, and non-zero, topological charge. 



2. g± from zero-mode wavefunctions 

In the £ -regime and in a fixed topological sector, correlation functions involving quark propa- 
gators may contain poles in 1/ (mV) n , where n is some integer number, whenever the contribution 
of the zero-modes to the spectral representation of the quark propagator gives a non-vanishing 
contribution to the correlation function. The residues of these poles are easier to compute than 
the correlation functions themselves. The idea, first explored in [||], is then to use the residues of 
the topological poles to perform the matching, instead of the full correlation function. Given a 
correlation function C v (x\,X2, ■■■), the residue can be isolated by 

Res 

Cy(x u x2,...) = , " +..., Res n = lim(mV)"C v (xi,x 2 ,...). (2.1) 

[my )" m— >0 

In ^ the two-point function of the pseudoscalar density was considered in this context. The 
presence of a pole in l/(mV) 2 implies that the corresponding residue can be computed fully in 
terms of the zero-mode wavefunctions, no propagator computation is required. On the effective 
theory side, the same pole does appear and the residue is a function of only the pseudoscalar 
decay constant, F, up to NLO. A numerical exploratory study in the quenched approximation was 
presented and the usefulness of the method to extract the low-energy coupling F was confirmed. 

In the present work, we extend this idea to the computation of three-point functions from 
which the weak low-energy couplings g± can be extracted. In particular we have considered the 
following ratios: 

Ra ^ lim m ^ (mV) 2 ^ y (d Xo P a (x)Q°(z)d yo P b (y)) v 

v ~ lim m _ K )(mV)i: x (o^P«(jc)/ig(z) ) v lim m ^ Q (mV)Zy{dy P b (y)J L b o(z)) v 
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Ay(x Q - zo,yp - go) + oA v (x - zq,yo - go) 
B v (*o- zo)fiv(>o -zo) 



cr = ±, 



(2.2) 



where P a = /*F7 5 r ax F, J L ^ = ^JqPiT"^ and 0* are the four fermion operators transforming in the 
84 and 20 representation of SU(4), for further details see If v,(jc) are the zero-mode wavefunc- 
tions of negative chirality, we find: 

A v (x Q -zo,y Q -z ) = 73 E (L nKx)S(x,z)rnP-Vi{z) £ ijJCy)5Cy,z)y M P_vj(i 



x.y.z igjf 



A v (xo-zo,Jo-zo) = ~ £( £ Tj)( JC )5( JC ,z)y M P_v ! -(z)T ] ;(y)%,z)y M P_v ; (z)) v; (2.3) 



x,y,z ijeJT 

where r\J(x)S(x,z) = l i^ x+a ^ s ^ x+a ^^^i^ x a ®) s ( x an( j s(x,y) is the quark propagator l . It is 
clear from eqs. ( |2~3| ) that a number of inversions equal to twice the topological charge, i.e. 2\v\, 
(with sources tj, (x) and T};(y), since xo and yo need to be fixed) is sufficient to construct the corre- 
lation function, whilst averaging over all the spatial positions of the three sources. Such averaging 
was only possible in the standard method of [Q] through low-mode averaging (LMA), and only for 
the contribution of the low-modes. The price of LMA is 12 + 2 x Ni uw inversions, where Ni ow was 
the number of low modes treated separately. Typically Ni ow can be as large as 20, and hence the 
numerical cost can be quite substantial. 



The matching of the amplitudes in eqs. ( |2.3| ) to the Chiral Effective Theory results in the 
following relation 



RGI 



z 2 



(2.4) 



RGI 



where fc ± (Miy) are the Wilson coefficients of the O operators, Z are their corresponding renor- 
malization factors and is the renormalization factor of the axial current. The RGI values of these 
factors have been computed non-perturbatively in [Jsj] . On the left-hand side of the equation, 3%^ is 



the ChPT prediction for the ratios of eq. (2.2) in the e regime. At NLO, only depend on F and 
the volume, as we show in the next section. 



3. Prediction in Chiral Perturbation Theory at NLO 



The result of the NLO computation of the two-point function in the denominator of eq. ( p.2| ) 
in the e expansion is 



(?) 



(*o) 



1 V 
l -l2 P -( 



;+*) 



FL) 



{FLf 



l*o 1 
T 



(3.1) 



where X =xq/T, p = T /L and is the number of dynamical flavours. The quenched result, 
is the same leaving out the term in parenthesis The first term corresponds to the LO result, 

that is the constant | v|. In Figure [I] we show the result for this quantity at NLO for different values 
of | v| in a symmetric box of size L = 2 fm. 



Similar expressions are obtained for the opposite chirality. 
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V=3 
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Figure 1: T3$ v (xq) for N — 2 (solid) and N = (dashed) as a function of xq/T. In both cases we take 
F =93 MeV and L = 2 fm. 



The results for the ratios of three to two-point functions are 



? v (xo-zo,yo-zo) 



V 



1±2 



F 2 ^V 



ply \gl(?x)+gl(?y)-hl(?x)-hi(T y )+lllTjj ) )H(T x ,T y ) 



where z x = (x - Zq)/T, i y = (y Q - Zq)/T, j8i is a constant [Q], and 



(3-2) 



iv i 

T| "2 "12 



«(t) = ww] 2 +i[A;w 2 +ipfA P (T) 2 ], 

//(T v , T y ) = ^ (T,)*i (%) - (T A . - T y ) - (/Jj (Z x ) - h\ {X x - Ty), 



(3-3) 

(3-4) 
(3.5) 



where \p\ = 2np\n\ with n = (n\ , «2 , «3 ) a vector of natural numbers, the primes refer to derivatives 
with respect to z, and 



cosh(|p|(|T|-i)) 
2|/?| sinh(|/?|/2) 



(3-6) 



Again the expressions in the quenched theory are identical at this order, except that the term (=F^) 
is not there. 

The ratios are therefore constants at LO, while at NLO they depend on the insertion of the 
sources. There is also a significant dependence on | v| already at LO. Both these features of the chi- 
ral corrections in differ from those of the ratios constructed out of left-current correlators, used 
in [Q], where no | v| dependence, nor temporal dependence was found at NLO. The large difference 
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between the chiral corrections in both cases probably implies different systematic uncertainties. 
Different systematics in the two approaches serve as additional consistency checks of the method. 

In the left plot of Figure | we show the result for M ±q v (T /3 - rT,2T/3 - rT)/(l T A) in a 
box of L = 2 fm for p = 1 and p = 2. Unfortunately NLO corrections seem to be rather large still 
at 2 fm, as shown on the right plot of Figure ||[ 

4. Exploratory quenched study 

We have recently carried out an exploratory study of the two and three-point functions in 
eq. (2^2) in the quenched approximation. We have considered the simulation parameters of Table [j], 



which correspond to two lattices with the same physical volume of around (2 fm) 4 , and different 
lattice spacings. 

We have used the overlap operator. For all the details on the implementation and algorithms 
we refer the reader to related previous work [{7|, ||, 0]. Here we will only present our preliminary 
results for the two-point functions. The results for the three-point functions will be presented in 
detail in a forthcoming publication. 

In the left plot of Figure ||, we show the results for B v (t) as a function of % = t/T compared 
with a fit of the form 

^ v (t) = a v + pJj-^\ , (4.1) 



where the NLO ChPT predictions of eq.(3.1) imply 



1 T I v 1 2 T I v 1 2 

av = |v| -T2lW ^ V = IW (4 " 2) 

We show in the right plot of Figure ||, the dependence of a v — \ v\ and j8 v on |v|, together with 
the fit to a parabola from which FL can be extracted. Although not very precise, the value of 
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Figure 3: Left: TB v (t) as a function of t/T for the lattice Al (full symbols) and Bl (empty symbols). 
The curves are fits of the Al data to NLO ChPT. Right: a v — \ v\ and j3 v as a function of |v| for lattice 
Al (full symbols) and Bl (empty symbols). The curves are fits to the NLO expectation for Al (solid) and 
Bl (dashed). 



Lattice 


P 


V 


|V| N conf 


x Q /a,y /a 


Al 


5.8458 


16 4 


1-5 282 


5,11 


Bl 


6.0375 


24 4 


2-5 236 


8,16 



Table 1: Simulation parameters 



F so obtained is in agreement with more direct determinations. Discretization effects in these 
observables are quite small. An interesting check of the two-point functions is provided by the 
chiral Ward identity, which implies the relation 

Z A B v (x -y ) = lim m 2 vY (P"(x) P°()0>v (4.3) 



The right-hand side of eq. (4.3) is the topological-pole contribution of the pseudoscalar two-point 
function, that is the observable that was studied in Since has been computed before 
we can compare the two sides of the equation. The ratio of the left over the right-hand side of the 
equation is shown in Figure || for the coarser lattice. The result is one in the chiral limit, as expected, 
in all topological sectors studied. Although it is no surprise that the Ward identity is satisfied 
given the exact chiral symmetry of the discretization, it is a non-trivial test of the method that it 
is also satisfied when both sides of the identity are truncated to the zero-mode pole contributions, 
especially because the limit m — ► can be done analytically on the right-hand side but not on the 
left-hand side. 
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Figure 4: Check of the truncated Ward Identity. 
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